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1 Introduction 



The notion of g-analogues (also called g-deformations) of soliton equations has been stud- 
ied from a variety of points of view. Kajiwara and Satsuma jl] obtained a g-analogue of 
the 1 + 1 dimensional Toda equation in a bilinear form. Partly motivated by this work, 
Mironov, Morozov and Vinet j2] introduced a g-analogue of the tau function of the Toda 
hierarchy. On the other hand, Wu, Zhang and Zheng j^] found a g-analogue of the KdV 
hierarchy along with soliton solutions. Frenkel |3] addressed this issue from the point 
of view of VT-algebras, and considered a g-analogue of the (generalized) KdV hierarchy 
and some other related soliton equaitons. Khesin, Lyubashenko and Roger (5j proposed 
a slightly different framework of g-pseudodifferential operators to formulate the Lax rep- 
resentation of these g-deformed soliton equations. Mas and Seco |E] applied the work of 
Khesin et al. to a class of g-analogues of W^-algebras. Adler, Horozov and van Moerbeke 
|l7|| pointed out that these g-analogues are closely related to the ordinary KP and Toda 
hierarchies. Having obtained a similar result in the context of bispectrality , Iliev con- 
structed a g-analogue of the KP hierarchy p. Tu pDj studied additional symmetries of 
Iliev's hierarchy. 

In this paper we consider a g-analogue of the modified KP (mKP for short) hierarchy. 
Our motivation is different from the preceding studies. Namely, we are interested in quasi- 
classical limit of the g-analogue as the underlying Planck constant h tends to 0. Actually, 
our true concern lies in the fate of the aforementioned g-analogue of the Toda hierarchy in 
the quasi-classical limit. The mKP hierarchy, originaly formulated as a system of bilinear 
equations connecting KP tau functions [TTl [T2], may be thought of as a subset of the 
Toda hierarchy |13jj. Because of this, one can use the mKP hierarchy as a prototype of 
technically more complicated consideration on the Toda hierarchy. 

In quasi-classical limit, the ordinary Toda hierarchy turns into the dispersionless Toda 
hierarchy ^3] • Accordingly, quasi-classical limit of the mKP hierarchy can be realized as 
a subsystem of the dispersionless Toda hierarchy. It is natural to expect a similar result 
for the g-analogue. To address this issue, we need a Lax formalism. We shall construct 
such a Lax formalism in the same framework as the ordinary Toda hierarchy. 

A few comments on the mKP hierarchy are in order. Firstly, various Lax representa- 
tions other than ours have been proposed for the mKP hierarchy ITTj ITHj. Secondly, 
the mKP hierarchy also appears in the work of Adler et al. [7] as the "one- Toda lattice" 
in their terminology. Thirdly, Takebe [I^ studied a generalization of the mKP hierarchy 
and its dispersionless limit. His construction is based on a different Lax formalism due to 
Dickey [H]. 
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This paper is organized as follows. Section 2 is a brief review on the tau function 
of the ordinary mKP hierarchy. In Section 3, its g-analogue is introduced and shown to 
satisfy a system of bilinear g-difference equations. In Section 4, these bilinear equations 
are translated to the language of wave functions, which turn out to satisfy a system of 
linear g-difference equations. These linear g-difference equations are used in Sections 5 
and 6 to formulate the Lax formalism and the description of quasi-classical limit. Section 
7 is devoted to the case of the g-analogue of the Toda hierarchy. 

2 Tau function of modified KP hierarchy 

The tau function r(s, t) of the mKP hierarchy depends on a discrete (Z valued) variable 
s and a set of continuous variables t = (ti,t2, . . .), and satisfies the bilinear equations 

/ r(s', t' - [X-'])t{s, t + [X-^])X''~'e^^'''''^^dX = (1) 

J A=oo 

for s' > s and arbitrary values of t' and t fUE]- The contour of the integral is understood 
to be a circle surrounding X = oo, and the following standard notations are used: 

oo 

[a] = (a, aV2, . . . , a'/k, . . .), ^t, A) = E ^.A^ 

n=l 

The bilinear equation for s' = s coincides with that of the KP hierarchy, so that T{s,t) 
is a tau function of the KP hierarchy. The bilinear equations for s' ^ s define a relation 
(a Backlund transformation in a generalized sense) connecting two solutions of the KP 
hierarchy. A pair of KP tau functions in such a relation is generated, for instance, by 
the action of a vertex operator jHlE]- If the chain of those KP tau functions t(s, t) is 
periodic, i.e., r(s + N,t) = r(s,t) for a positive integer A^, the relation reduces to the 
ordinary cyclic Darboux transformations for the A^-th generalized KdV hierarchy. In the 
simplest case of A^ = 2, the reduced hierarchy with two tau functions r(0, t) and r(l, t) is 
nothing but the modified KdV hierarchy. This explains the origin of the word "modified" 
in the name of the hierarchy. 

Actually, yet another "modified" set of bilinear equations can be derived from these 
fundamental bilinear equations by shifting t' as 

t'^t'-[X^'] [X7A 

Xi, . . . ,Xn being assumed to be on the far side (i.e., in the domain that contains A = oo) of 
the integration contour. The exponential function e^*^*'~*'^^ thereby gives rise to an extra 
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factor 

so that the bihnear equations take the "second-modified" form 

/ T{s',t'-[Xi'] [x--]-[x-'Ms,t+[x-']) 

J X=oo 

X A^'-^e«(*''*'^)(A - Ai) ■ ■ ■ (A - Xn)dX = 0. (2) 

If Aj's are suitably chosen, as we shall see, these bilinear equations turn into g-difference 
equations. 

3 Q'-analogue of mKP tau functions 

We now introduce a new set of continuous variables x = {xi, X2, ■ ■ ■) and a set of parame- 
ters q = (gi, g2, • • •) with |gi| < 1, \q2\ < 1, etc., and consider the multivariate deformation 



T,{s,t,x)=Tis,t+Y.[Xnt^) (3) 
n=l 



as a g- analogue of the foregoing mKP tau function. In this definition, [a]^"'' stands for a 
g-analogue of the [a] symbol of the form 

[«r-fo,...,0,a,0,...,0,%Zi)!^,...,0,..^^ ^ 



2(1 - g2) A;(l - g^ 



with the non-zero elements (1 — q)^a^ /k{l — g^) placed in the fcri-th component for k = 
1,2,.... This is essentially the same change of variables as proposed by Mironov et al. j2] 
for the case of the Toda hierarchy, except that we keep the old variables ti, t2, ... along 
with the new ones xi, X2, .... If the XnS other than xi are set to zero, Tq{s,t,x) reduces 
to the univariate deformation 

T,{t,x)=T{t+[x],), [X],= [X]^^\ (4) 

of the KP tau function r(t) studied by Adler et al. [7' and Iliev [H] in their work on a 
g-analogue of the KP hierarchy. 

is defined so as to satisfy the g-difference relation 

[HS"^ = [«r-[(i-9)«]^"\ (5) 
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where [a] '■"'^ denotes a modification of [a] of the form 

f0,...,0,a,0,...,0,y,...,0,...,0, 

in which the fc-th entry jk of [a] is relocated to the fcn-th component and the other 
components are set to zero. This is a generahzation of the g-difference relation of 
that lies in the heart of the work of Adler et al. [Zj and Iliev . As noted by Mironov et 
al. [a]*^"-* can be expressed as a linear combination of [a] symbols, 

n 



so that the g-difference relation of [aj^"^ can be rewritten as 

n 

Mi"^ = Mg"^ - E[e'"'^"(l - (6) 
i=i 

This g-difference relation enables us to relate the multiplication Xn —>■ qnXn in rq(s, t, x) 
to the shift t^t - [A^^] [A;^] in ©. If A/s are indeed set to the values 

A^. = e-2-^J'/"(l - g„)"'/"a;;^/" 

that can be read off from the last polynomial factor in Q takes the remarkably simple 
form 

(A - Al) ■ ■ ■ (A - A„) = A" - (1 - qnr'x-\ 
Thus the modified bilinear equation (0) can be converted to the bilinear equation 

Tq{s\t' - [A^^], . . .,qnXn, ■ ■ .)Tg{s,t+ [A"^], . . . , . . .) 

1 

X A^'-^e«(*'-*'^) (1 - (1 - g„)x„A")c/A = (7) 
for Tq{s,t,x). Moreover, iterating © k times gives higher order g-difference relations 



A=oo 



n k—1 



from which a similar bilinear equation can be derived. Actually, one can consider simul- 
taneous g-multiplication of (but fc„ = except for a finite number of 
n's), which leads to bilinear equations of the more general form 



/ r^is',t'-[X-\x')T,{s,t+[\-%x)\ 

J A=oo 



OO kn~l 

X n 11(1 - (1 - qr,)qtxnXndX = 0, (9) 



where x' = {x[, x'a, . . .) and x = {xi,X2, ■ ■ ■) are related by nonnegative powers 92^, • • • 
of qi,q2, ... as 

x[ = g^xi, 4 = q^'x2, .... (10) 



4 (/-analogue of mKP Wave functions 

Since the time variables t of the mKP hierarchy are retained in the definition of Tq{x, t, x), 
one can define an associated wave function \l/g(s, t, x, A) and its conjugate t, x, A) as 
wave functions of the mKP hierarchy: 

Tq[S, t, X) 

n',t,x,X) = ^''^''';^[^'y\ -^e-^^'''^e,{x,Xr\ (H) 
Note that the exponential part A^^e^^*^*'^) is multiplied by the extra factor 

oo 

e,{x,X) = IleZ""^ (12) 

n=l 

where denotes the so called g-exponential function 

< = exp E un 1 = n (1 - (1 - ^W)-'- 



k=0 



This extra factor originates in the shift of t by the sum of [xn]^"^ in the definition of 
Tq{s,t,x), and behaves exactly like the ordinary exponential factor under the action of 
g-difference operators. 

A prototype of this g-exponential factor is the g-exponential function e^^ that plays a 
central role in the g-analogue of the KP hierarchy El . It satisfies the g-difference 
equation 

Dq{x)ef = Aef (13) 
for the g-difference operator Dq{x) that acts on a function of x as 

Dq{x)f{x) = . 

qx — X 

Let us also mention that the usual Leibniz rule is modified to this operator as 

Dq{x){fg) = Dq{x)f ■ g + Tq{x)f ■ Dq{x)g, (14) 
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where Tq{x) denotes the g-multiphcation operator 

Tg{x)f{x) = f{qx). 

By the same token, as imphcitly noted in the work of Mironov et aL the q- 
exponential factor eq{x, A) satisfies the g-difference equation 

Dq^{xn)eg{x, A) = A"eq(x, A) (15) 

or, equivalently, 

Tq„ixn)eg{x, A) = (1 - (1 - g„)x„A")eg(x, A). (16) 
Iterating the last equation yields 

Tg„{xn)''"eqix, A) = n (1 - (1 - Qn)qi"XnX'') ■ eq{x, A). 



The prefactor on the right hand side coincides with the last polynomial factor in (jH)). This 
means that (jH)) is actually a bilinear equation for the wave functions: 

(f ^,(s',t',x',A)^;(s,t,x,A)rfA = 0. (17) 

J A=oo 

We can now follow a standard procedure to derive a system of linear g-difference 
equations for "^q^s, t, x, A) from the bilinear equation (fTTjl . A technical clue is the following 
lemma. 

Lemma 1 If ^{s,t,x, X) is a function (or Laurent series) of the form 

oo 

$(s, t, x,\) = J2 Ms, t, x)A-" ■ A'e«(*'^)eg(x, A) (18) 

n=l 

and satisfies the bilinear equation 

(f ^{s',t,x,X)'^q{s,t,x,X)dX = (19) 

J A=oo 

for s' > s and arbitrary values of (t, x), then $(s, t,x, X) = 0. 

Proof. The bilinear equation in the statement of the lemma implies the equations 

s'-s+l 

J2 (Pn{s', t, x)m;*,_,_„+i(s, t,x) = 

n=l 



for (f)n{s,t,x)'s and the coefficients of the Laurent expansion 

oo 

%{s,t,x,X) = (l + 5:<(s,t,x)A-")A-^e-«(*'^)e,(x,A)-^ 

n=l 

of \I/*(s, t, X, A). Starting with the case of s' = s where this equation reduces to 0i(s, t, x) = 
0, one can easily show by induction that (f)n{s, t,x) = for all n. 

To derive a linear g-difference equation for ^^(s, A), let us note that the action 
of Dg^{xn) on "^gis, t, X, A) yields a function of the form (A" + ■ ■ ■)ye^^^'^^eg{x, A), where 
A + - ■ ■ is a Laurent series of A that starts from A" and continues to terms of A*^"^, A"~^, . . .. 
One can find a difference operator Cn (i.e., a linear combination of the shift operators 
^mds -j-j-^g^^ g^^^ ^ function of s as e™'^=/(s) = /(s + m)) of the form 

n 

Cn = e"^= + J2 Cnm(s, t, x)e("-"^)^= (20) 

m=l 

that absorbs the main part of the Laurent expansion of Dg„(x„)\E'g(s, t, x, A) as 

Dg^ {xn)^g{s, t, X, A) - Cn^gis, t, X, A) = 0(A-i)A^e«(*'")e,(x, A). 
Let $(s,t,x, A) denote this quantity. One can derive, from (|T7jl . such relations as 
/ Tg„(x„)^'g(s,t,x,A) ■ ^l{s,x,t,X)dX = 

J A=oo 

and 

<f e"'^^^g{s,t,x,X) ■ %{s,t,x,X)dX = 

J A=oo 

for m > 0. These relations imply that the assumption of the aforementioned lemma holds 
for ^{s,t,x, X), which thereby turns out to vanishes identically. One can thus conffim 
that \E'q(s, t, X, A) satisfies the linear g-difference equation 

^g„(a;„)*g(s,t,x,A) = C„^g(s, t, x. A). (21) 

5 g'-analogue of Lax, Zakharov-Shabat, and Sato equa- 
tions 

Consistency of the linear g-difference equations ()21|) leads to a system of zero-curvature 
equations for C„'s. To derive the consistency condition, we start from the commutativity 
relation 

Dq„ {Xn)Dg^{Xm)'^g{s, t, X, A) = (x„)Dg„(x„)*g(s, t, X, A), 
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and calculate both hand sides using (PT|) and the g-difference Leibniz rule. The outcome 
is the equation 

which implies that the g-difference analogue 

-Dq„{Xn)Cm " Dq^(^Xm)Cn + Tq^(Xn)Cm ' Cn ~ Tq^^^(^Xm)Cn ' Cm = (22) 

of the zero-curvature equation (also called "Zakharov-Shabat equations" ) is satisfied. 
We now introduce the dressing operator (also called "Sato- Wilson operator" ) 



W = l + Y.Wnis,t,x)e-'"'% (23) 

n=l 

where w„(s,t,x) denote the coefficients of the Laurent expansion 

oo 

^qis, t, a;, A) = (l + ^ w^is, t, x)\-^) A^e«(*'^)e,(x, A). 

n=l 

(PT|) can be converted to the g-difference equations 

D,MW = CnW - Tq„{xn)W ■ e"^^ (24) 

One can eliminate Cn to derive a nonlinear system of equations for W. This is achieved 
by multiplying both hand sides by W from the right side, 

D,MW ■ = Cn- TqSxn)W ■ e^^^W-\ 

and separating the difference operators into the part of nonnegative/negative powers of 
e^" . Let (A)>o and (A)<o denote the projection of a difference operator to these parts: 

n — n>0 n n<0 

The (■)>o part of the foregoing g-difference equation gives the relation 

Cn={TqSxn)W-e^^^'W-')^^. (25) 

This enables one to eliminate C„ from the g-difference equation itself. The outcome is 
the nonlinear g-difference equation 

D,Sxn)W = -{TqSxn)W ■ c'^'^W-') ^W. (26) 
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This is a g-difference analogue of the so called "Sato equations" . 
Finally we introduce the Lax operator 



L = We^^W-' = e^^ + J2 Une^'~"^^% (27) 

n=l 

which turns out to satisfy the g-difference Lax equations 

Dg^iXn)L = CnL - Tg„(x„)L ■ Cn (28) 



as a consequence of the Sato equations (j26|) . The same Lax equations can also be derived 
from the consistency of the formal spectral equation 

L^g(s,t,x, A) = A^,(s,t,a;, A) (29) 

with the linear g-difference equations (PT|) . The commutativity, Dq^^{xn)Dq^{xm)L = 
Dq,^{xm)L, of the flows of the g-difference Lax equations is ensured by the g-difference 
Zakharov-Shabat equations (j^^ . 

To conclude, let us mention that t, x, A) satisfies the linear equations 

9t„^,(s,t,x,A) = 5„^g(s,t,x,A) (30) 

of the ordinary mKP hierarchy as well, -B„ being a difference operator of the form 

Bn = e"^^ + J2 bnmis, t, x)e("-"^)^= . (31) 

m=l 

One can derive from these linear equations the Sato equations 

dt„W = -{We''^'W-^)^W, (32) 
the Zakharov-Shabat equations 

dt„Bm — dt^Bn + [Bm, = 0, (33) 



and the Lax equations 



along with the relation 



dt„L=[Bn.L\ (34) 



B^=(We^^^W-^)^^ = {L^)>o (35) 

connecting Bn with W and L. These equations comprise a Lax formalism of the mKP 
hierarchy. What we have derived above is a g-difference analogue thereof. 
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6 Quasi-classical limit and Hamilton-Jacobi equations 



We now turn to the issue of quasi-classical limit. To this end, we set the parameters 
to depend on the Planck constant h as 

q^ = q'' = e-""^^, q = e-^^ (36) 

where (3 is an arbitrary constant, and rescale the independent variables s,tn,Xn as 

S ^ s/h, t„ tn/h, Xn XnJ (1 - (37) 

The fundamental operators dt,^ and D thereby rescaled as 

which resembles the usual set-up of quantum mechanics. The linear equations of the wave 
functions are reformulated in the ^-dependent form as 

hdt„^,{X) = 5„^,(A) (38) 

and 

(l-g")D,n(a;„)^,(A) = C„^,(A). (39) 

Also note that we have slightly changed the notation, namely, we omit writing the de- 
pendence on X, t, X, h to highlight the dependence on A. The coefficients bnm and c„m 
of 



TD _ nhds I \^ U Jn-m)nds fi _ nhds i ST^ „ pi 



^{n—m)hds 

m=l m=l 



too, are understood to depend on h arbitrarily. To achieve reasonable quasi-classical limit, 
however, we assume that bnm and c„m have a smooth limit 

lini6„^ = 6°^, lini c„^ = c°„ (40) 

as h 0. As in the case of the Toda hierarchy |l4j, this imposes a rather strong condition 
on the solution in consideration; let us simply assume that these conditions are satisfied. 
One can readily see from ()25|1 and that i?„ and C„ coincide in the limit as h —>■ 0, 
namely, 

6° = c° f411 
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Under these assumptions, one can set the wave function in the WKB form 

^^(A) = exp{h-^SiX) + 0(/i°)) (42) 

or, more loosely. 

Since the g-exponential factor eq{x, A) behaves as 

e,(x. A) = exp (^f ^^Mj^ + 0(^°)) , (43) 
where Li2(-2) denotes the dilogarithmic function 

oo 

Li2(2:) = E p ' 
k=i 

the phase function S{X) = S{s,t,x, A) is a Laurent series of the form 

S{X) = s log A + e(t. A) + ^ " ^ + (negative powers of A). (44) 

n=l ''^P 

One can derive from ()38j) and a system of Hamilton- Jacobi equations as follows. 

For warm-up, let us first consider (p?Hj) . Actually, this has been done in the study of 
quasi-classical limit of the Toda hierarchy |14 . Upon substituting the WKB ansatz, both 



hand sides of can be expressed as 

dt^^,{X) ~ 9t„5(A)e^"^W, B^^,{X) ~ S„(e^=^(^))e^"^W, 
where Bn{p) is a polynomial of the form 

n 

^n(p)=p"+EO"-™ (45) 
m=l 

Thus one obtains the Hamilton- Jacobi equation 

d,,^SiX) = i3„(e^^^W). (46) 

Let us now consider ()39|) . The main task is to determine the asymptotic form of the 
left hand side, i.e., 

(l-g")D,Jg")vI/,(A) 
= x-'[exp{h-'S{X) + 0{h')) - exp{h-'T,,Xq'')S{X) + o{n'))). 
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Ti.,^(g")S'(A) can be expanded to a Taylor series of Ti as 

T,S(f)S{s,t,x,\) = 5(s,t,...,e-"^^x„,...,A) 

= S{s, t, X, A) — n(3hxndx„S{s, t, x, A) + 0{h^). 

Consequently, 

(1 - g")Z}.„(g")vl/,(A) ~ x-\l - exp(-n/5x„9.„5(A)))e^"'^W. 
On the other hand, since Cnm and hnm coincide in the limit as — > 0, 

Thus one obtains the Hamilton- Jacobi equation 

x-\l - exp(-n/3x„9,.„5(A))) = K„(e^^^(^)). (47) 
It will be also suggestive to rewrite this equation as 

Xnd.„S{\) = log(l - x„S„(e^=^("))). (48) 

One can now follow the case of the Toda hierarchy [14^ to convert these Hamilton- 
Jacobi equations to a Lax formalism of the "dispersionless" type ^H]- The Lax function 
C{j)) is obtained by solving the equation 

dsS{\)=p (49) 

for A as 

oo 

\ = c{p) = p+Y.uy-\ (50) 

n=l 

The coefficients vP^ coincide with the limit of the coefficients m„ of L as ^ ^ 0. ()46p and 
(I48|l correspond to the Lax equations 

dt^p) = {Bn{p),C{p)}, 

Xndx„C{p) = -{\0g{l -XnBn{p)),C{p)}, (51) 

with respect to the Poisson bracket 

{/, 9} = pdpf ■ ds9 - dsf ■ pdpQ. 
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One can further introduce the Orlov-Schulman function 



M{p) = XdxS{X)\ X=C(p), (52) 
which turns out to satisfy the Lax equations 

dt^Mip) = {Br.{p),M{p)}, 

Xnd.„M{p) = ~^{\og{l -xMp)),M{p)} (53) 
np 

and the Poisson commutation relation 

{C{p),M{p)} = C{p). (54) 

7 g-analogue of Toda hierarchy 

The tau function T{s,t,t) of the Toda hierarchy [T3] depends on yet another set of con- 
tinuous variables i = (ti,t2, • • •)' ^^"^ satisfy the bilinear equations 



A=oo 



t{s', t' - [\-\i')T{s, t + t)A^'-^e«(*'-*'^)dA 



A=0 



-{s\t\t' - [A])r(s,t,t + [A])A^'-^e«(*'"*"^"')dA, (55) 



where s' and s ar now aribrary (namely, there is no ordering constraint), and the contours 
of the integrals on both hand sides are understood to be a circle surrounding \ = oo 
and A = respectively. Following Mironov, Morozov and Vinet 0, we now introduce 
two sets of continuous variables x = {xi,X2, . . .), x = {xi,X2, ■ ■ ■) along with parameters 
q = (gi, g2, • • •)) ^ = (^i; ^2, • • •)> consider the g-analogue 

oo oo 

Tg,g{s,t,t,X,x) = t{s, t + E Hi?' ^+ J2[^n& (56) 

n=l n=l 

of r(s, t,t). The modified bilinear equations Q can be generalized to Tq^g{s, t, i, x, x). One 
can rewrite those modified bilinear equations for Tq^q{s, t, i, x, x) to the bilinear equations 

/ ^g,5(s', t', i', x', x', A)^;g-(s, t, t, X, X, A)dA 

J A=oo 

^q,q{s', t', F, x', x', \)^* g{s, t, t, X, X, \)d\ (57) 



/A=0 

for the wave functions 

^q,q{s,t,i,x,x,X) = ^^•^"^'''~f-'^'^''^'^^ -e^(^'^)e,(x,A), 
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A) = ^5 - 1,M+ [^U- ^) A-^e-^(U'%,(,;, A-^)-^ (58) 

x'j^,x'^,Xn,Xn are understood to be related by nonnegative powers g^" and g^" of g„ and 
Qn as 

x'l = q'l^xi, X2 = qi'x2, x[ = g^Xi, 4 = g^^^a, ■ ■ ■ ■ (59) 

The wave functions turn out to satisfy a set of linear differential and g-difference 
equations. To this end, the previous lemma has to be generalized to an equation of the 
form 

^(s', t, t, X, X, X)^q^q{s, t, t, X, X, X)dX 
J \=oo 

= i ^{s',t,t,x,x,X)^g^g{s,t,t,x,x,X)dX, (60) 

where $(s, t, t, x, x, A) is the same as in the previous case and <l>(s, t, t, x, x, A) is a function 
of the form 

oo 

^(s, t, t, X, x,X)^Yl Ms, t, t, X, x)A" • A"e^(*'^"')eg-(x, A"^). (61) 

n=l 

The generalized statement is that such a pair of functions $(s, t, i, x, i, A) and $(s, t, i, x, i, A) 
vanish identically. One can threby show that ^^,5(5, t, i, x, x, A) satisfy the linear q- 
difference equations 

Dg^{Xn)'^q,q{s,t,t,X,X,X) = Cn"^ g,q{s , t,t, X , X, X) , 

Dqr,{Xn)^q,q{s,t,t,X,X,X) = C^^ q^q{s, t,t, X , X, X) (62) 

along with the linear differential equations 

dt„'^q,q(s,t,t,X,X,X) = Bn'ifq,q(s,t,t,X,X,X), 

dtn'^q,q{s,t,t,X,X,X) = Bri'^q^q{s,t,t,X,X,X) (63) 

of the ordinary Toda hierarchy, and that the same equations hold for ^g^g{s,t,t,x,x, X) 
as well. Bn and C„ are the same difference operators as in the previously case. B^ and 
Cn are difference operators of the form 

n— 1 n—l 
m=0 m=0 
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Quasi-classical limit is achieved by setting qn = Qn = q"", Q = e~^^, and rescaling 
the independent variables as s ^ s/h, tn tn/h, t„ ^n/^, a^n Xn/{1 — q""), Xn — > 
Xnl (1 — q^)- Also the coefficients of i?„, C„, Cn are allowed to depend on Ti arbitrarily, 
except that they have a smooth limit as ^ 0. Assuming the WKB ansatz 

^,,(s,A) ~exp(;i-iS(A) + 0(/i°)), (65) 

one can derive the Hamilton- Jacobi equations 

Xnd^^SiX) = -i-log(l-S„(e^=^W)), 
np 

Xnd.^SiX) = -J-log(l-g„(e^=^W)), (66) 

np 

where Bn{p) is the same polynomial as in the case of the modified KP hierarchy and Bn{p) 
is a polynomial in of the form 

n-l _ 

^n(p) = E O'"-", 6°^ = lim &_. (67) 

m=0 

One can define the Lax function C{p) and the Orlov-Schulman function Ai{p) in the same 
way as the previous case. They satisfy the Poisson commutation relation 

{C{p),M{p)} = C{p) (68) 

and a system of Lax equations with respect to the same Poisson bracket {-,•}• The Lax 
equations for C{p) read 

dt^p) = {Bn{p),C{p)}, 
d,„C{p) = {B^ip),Cip)}, 

Xndx„C{p) = - ^{log{l - XnBn{p)),C{p)}, 

Xndx„C{p) = ^ {log{l - XnBn{p)),C{p)}, (69) 

np 

and Lax equations of the same form hold for Ai{p). One can repeat the same calculations 
for the WKB ansatz 

^,,,-(A) ~ exp{h-'S{X) + 0(/i°)) (70) 

of 'l'g^g(A). This eventually leads to the second pair C{p), Ai{p) of Lax and Orlov-Schulman 
functions. Using the four functions C{p), Ai{p), jC{p), Ai{p), one can formulate a "twistor 
theory" [TB] of the extended Lax formalism. 
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8 Conclusion 

We have considered a g-analogue of the mKP hierarchy as a prototype of the g-analogue 
of the Toda hierarchy. As expected, we have been able to describe the quasi-classical 
limit of both systems by the method developed for the Toda hierarchy 14J. An unusual 
feature is that the g-difference flows turn into the flows with logarithmic generators such 
as log(l — XnBn{p)) and log(l — This is reminiscent of a logarithmic structure 

that can be seen in Takebe's result [THj . 

As a final remark, let us mention a possible link of the present work with the random 
partition calculus in gauge theories and topological strings pUl I21j . Such a link can be 
seen in the simple (almost trivial) tau function 



of the Toda hierarchy, where A is an arbitrary constant. If g„, g„ are set to the special value 
q"', and all tn,in,Xn,Xn but xi and Xi are turned off, then the g-deformed tau function 
reduces to 



which, upon suitably adjusting the variables Xi,Xi, coincides with a relevant physical 
quantity (such as the partition function) in the random partition calculus. 
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